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Abstract. Using the intertwining relation we construct a pseudosuperpartner for a (non-Hermitian) Dirac- 
like Hamiltonian describing a two-level system interacting in the rotating wave approximation with the electric 
component of an electromagnetic field. The two pseudosuperpartners and pseudosupersymmetry generators close 
a quadratic pseudosuperalgebra. A class of time dependent electric fields for which the equation of motion for a 
two level system placed in this field can be solved exactly is obtained. New interesting phenomenon is observed. 
There exists such a time-dependent detuning of the field frequency from the resonance value that the probability 
to populate the excited level ceases to oscillate and becomes a monotonically growing function of time tending 
to 3/4. It is shown that near this fixed excitation regime the probability exhibits two kinds of oscillations. The 
oscillations with a small amplitude and a frequency close to the Rabi frequency (fast oscillations) take place at 
the background of the ones with a big amplitude and a small frequency (slow oscillations). During the period 
of slow oscillations the minimal value of the probability to populate the excited level may exceed 1 /2 suggesting 
for an ensemble of such two-level atoms the possibility to acquire the inverse population and exhibit lasing 
properties. 
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1. Introduction 

The supersymmetry in physics has been introduced in the quantum field theory for unifying 
different interactions in a unique construct pQ. Supersymmetric formulation of quantum 
mechanics is due to the problem of spontaneous supersymmetry breaking [2\. Ideas of 
supersymmetry have been profitably applied to many nonrelativistic quantum mechanical 
problems since, and now there are no doubts that the supersymmetric quantum mechanics 
(SUSY QM) has its own right to exist (for recent developments see a special issue of Journal of 
Physics A, vol. 34, No 43, 2004). It is worth noticing that most papers in this field deal with 
the Hermitian Hamiltonians. 

Differential equation of Schrodinger-like type with a non-Hermitian Hamiltonian appears 
in many physical models. One can cite quantum systems coupled to the environment like a 
hydrogen "atom" in an interacting medium subject to a dissipative force |3] (see also or 
different decay or collision reactions (see e. g. j^j; for more recent developments see 0; in [7] 
the method of SUSY QM is involved). Physical needs initiated a deep mathematical study of 
spectral problems with non-Hermitian Hamiltonians in 50th and 60th of the previous century. 
The most essential result was first obtained by Keldysh [8 who proved the completeness of 
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the set of eigenf unctions and associated functions for a regular Sturm-Liouville problem with 
a non-Hermitian Hamiltonian. In the books by Naimark [5] and Marchenko ^U] one can find 
good reviews of these studies. 

A new impact to studying different properties of non-Hermitian Hamiltonians is due to 
the discovery that the real character of the spectrum of a non-Hermitian Hamiltonian may 
be in particular related with so-called PT-symmetry jTT] and suggestion to generalize the 
quantum mechanics by accepting non-Hermitian Hamiltonians with a real spectrum to describe 
physical observables ^2] (see also the review [T%]). The necessary condition for such a 
generalization consists in the possibility to define a Hilbert space with a positive definite 
metric which is intimately related with the property of a Hamiltonian to be diagonalizable 
(for recent discussions see e.g. [HJ ITH]). This apparently may be assured in many cases 
since non-diagonalizable Hamiltonians may be transformed into diagonalizable ones by SUSY 
transformations JO] • The latter property permits us to suppose that the method of SUSY QM 
may become an essential ingredient of the complex quantum mechanics. This conjecture is also 
supported by established properties of this method not only to offer possibility for obtaining new 
exactly solvable complex potentials from known ones but also to help deeper understanding 
different properties of complex potentials [13 US]. In particular, an explicit construction of a 
superalgebra involving non-Hermitian Hamiltonians, which may be useful in different contexts 
i.e. integrability, quantization, different quantum field models etc, is shown to be possible 
[TU] and even is now developed till the notion of pseudosupersymmetry [2Q. and nonlinear 
pseudosupersymmetry [2T] . 

The relation of the general two-level model described by a non-Hermitian Hamiltonian 
acting in the two-dimensional Hilbert space C 2 with the pseudosupersymmetry is discussed 
by Mostafazadeh [20] . In contrast to the approach of this author we reduce the time-dependent 
Schrodinger equation for the two level system, interacting in the rotating wave approximation 
with the electric component of an electromagnetic field, with a Hermitian Hamiltonian (see 
e.g. [22]) to the one- dimensional stationary Dirac equation with an effective non-Hermitian 
Hamiltonian where the time plays the role of the space variable. If we considered the spectral 
properties of the latter Hamiltonian we would define it in the Hilbert space L 2 (0, T) <g> C 2 . But 
as we shall see in our approach the spectral parameter in the Dirac equation is not related 
with spectral properties of the two-level system. Therefore we will not discuss any spectral 
features of this Hamiltonian and in particular its diagonalizability. Of course, the obtained 
Dirac equation is completely equivalent to the initial Schrodinger equation and if one studied it 
by usual means one would not get any new information about the two-level system. From this 
point of view the method of SUSY QM we are using proves its extreme efficiency once again. 

To find a pseudosuperpartner for the given Dirac-like Hamiltonian we are using the technique 
of intertwining operators developed in [23] for the one- dimensional stationary Dirac equation. 
We have to notice that the application of results of this paper to our particular problem is 
not straightforward since transformation operators of the general form do not preserve the 
very peculiar form of the effective Dirac Hamiltonian corresponding to the two-level system. 
So, below we show how from the wide variety of possible transformations one can choose 
the necessary ones. In our approach in contrast to [20] the two pseudosuperpartners and 
pseudosupersymmetry generators constructed with the help of first order intertwiners close a 
quadratic pseudosuperalgebra. As it usually happens for the method of intertwining operators 
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[21] if one of the two Hamiltonians is exactly solvable the same property takes place for the 
other. In this way starting from the simplest case corresponding to the famous Rabi oscillations 
we have found new electric fields having time-dependent frequencies for which the equation of 
motion of the two-level system has exact solutions. While analyzing solutions of the Schrodinger 
equation we have found a new interesting physical phenomenon. We show that there exists such 
a time-dependent detuning of the field frequency from the resonance value that the probability 
to populate the excited level ceases to oscillate and becomes a monotonically growing function of 
time tending to 3/4. Of course this is a strictly fixed excitation regime similar to resonance. We 
also study how the above probability behaves under small deviations from this specific regime. 
We have found that when the parameters of the model are close enough to the specific values 
the probability exhibits two kinds of oscillations. The oscillations with a small amplitude and 
a frequency close to the Rabi frequency (fast oscillations) take place at the background of the 
ones with a big amplitude and a small frequency (slow oscillations). During the period of slow 
oscillations, which grows when the parameters of the model approach the above specific values, 
the minimal value of the probability to populate the excited level may exceed 1/2 suggesting 
for an ensemble of such two-level atoms the possibility to acquire the inverse population and 
exhibit lasing properties. 

We have to notice that some of the results we expose below are known from the previous 
paper [23] • These authors also use a similar intertwining technique but they do not relate 
it with the pseudosupersymmetry and do not give any analysis of solutions this method can 
provide with. Moreover, we give a deeper analysis of restrictions imposed on transformation 
operators by the features of the two-level system. In particular, we show that both the new 
Hamiltonian and solutions of the new Dirac equation can be expressed in terms of a real- valued 
function which is a solution of a second order differential equation with real coefficients. Since 
such equations have real solutions always our analysis opens the direct possibility to realize 
chains of transformations preserving the form of the Dirac-like Hamiltonian imposed by the 
features of the two-level system. 

2. Preliminary 

The two-level model in the rotating wave approximation with a possibly time-dependent 
detuning is described by the following system of equations (see e.g. [22]): 

iA 1 - fA 1 = £A 2 iA 2 + fA 2 = £Ai (1) 

where £ = ^E di 2 , d\ 2 is the matrix element of the dipole interaction operator, E is 
the amplitude of the electric component of an external electromagnetic field; / = h4?{8t), 
S(t) = 012 — 0(i), 012 = \{ £ i ~ £ 2), £ i and e 2 are energy levels of the free atom and 0(t) is 
the field frequency; the dot over the symbol means the derivative with respect to time. While 
normalized properly the functions | ^4.! (t ) | 2 and ^(t)] 2 give occupation probabilities for the 
ground and excited states respectively. If does not depend on time (hence f = \5 = const) 
solutions of the system ([!} are well-known. For instance, with the initial condition A 2 = and 
A\ = 1 at t = we get the well-known formula [22J for the excited state occupation probability 
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if initially the system is in the ground state 

Pit) = \A 2 (t)\ 2 = ^5 [1 - cos(20i)] fi 2 = f 2 + e (2) 

with 2£ known as the Rabi frequency. The probability (j2J) is an oscillating function of time (so 
called Rabi oscillations). At the resonance (/ = \8 = 0) it oscillates with the Rabi frequency. 
Therefore the value 5(t) characterizes the detuning of 0(t) from its resonance value equal 0x%. 
In Section 5 using the formalism developed in Section 4 we shall get time-dependent functions 
/ = f(t) (and hence S(t)) for which system (0) permits exact solutions. As we show below 
(Section 5) time-dependent corrections to the detuning that we will consider although may 
change crucially the time-dependent behavior of the solutions of system but they essentially 
keep oscillating character of the probability to populate the excited level with the frequency 
close to 2Q. Yet, the absence of the Rabi oscillations may be considered as oscillations with 
the same frequency but with the zero amplitude since they may be obtained as corresponding 
limiting case of oscillations with a non-zero amplitude. So, in our approach the rotating wave 
approximation is as good as it is in the classical case of the electric field of a constant frequency. 

Let us rewrite system (|1J) in the matrix form 



where 



h q = h = ryd t + V (3) 



V = ifoay (4) 



7 = ia x , E = £, ^ = (A 1 ,A 2 ) T (the superscript "T" denotes the transposition) and we 
replaced / (which we will call the "potential") in (JTJ) by / ; cr x ,y,z denote the standard Pauli 
matrices. Equation (jSJ) is the one- dimensional stationary Dirac equation with the non-Hermitian 
Hamiltonian h Q defined by the potential (J3J) where t plays the role of the space variable. By 
construction the parameters /o and E are real. For a fixed value of the dipole momentum of 
the irradiated system the parameter E = £ is defined by the amplitude of the electric field and, 
hence, is not related with spectral properties of the system. A useful comment is that since the 

( f £ \ 

Hamiltonian of the system (0) is Hermitian, H sc h — I c f ) , the evolution of the two-level 



system is unitary even for a time-dependent function / = f(t). This means that the C 2 inner 
product, |Ai(t)| 2 + |y4 2 (t)| 2 , for the Dirac equation (jSJ) is t- independent. 



3. SUSY algebra with non-Hermitian Hamiltonians 



Let us have a non-Hermitian Hamiltonian ho. We will not consider it as a Hamiltonian acting 
in a Hilbert space but to construct a SUSY algebra we need adjoint operators which we will 
introduce in a formal way. Denote by Hq the operator formally adjoint to ho- As usual the 
adjoint operation consists in taking the complex conjugation and transposition, the operator 
of the first derivative is skew- Hermitian and (AB) + = E> + A + . 

Let h% be a "transformed Hamiltonian" which should be found together with the 
transformation operator L by solving the intertwining relation Lh = h\L and h\ be its adjoint. 
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The later participates in the adjoint intertwining relation h$L + = L + hf. It means that the 
operator L + transforms eigenfunctions of H{ into eigenfunctions of Hq . 

Let us suppose that there exists an operator J such that h^ x = Jh 01 J and J 2 = ±1, 
J + = ± J (in general both signs may be accepted). Then from the adjoint intertwining relation 
it follows that JL + Jh\ = h$JL + J meaning that the operator JL + J realizes the backward 
transformation from h\ to ho and the operator JLJ transforms from Hq to h x . From here we 
infer that the superposition JL + JL transforms solutions of the equation (J3J) into solutions of 
the same equation meaning that this is a symmetry operator for this equation. In the simplest 
case when L is a differential operator that we would like to consider this symmetry operator 
may be a function of h , so we will suppose that JL + JL = Fi(h ). By the same reason the 
superposition LJL + J may be a function of hi leading to LJL + J = F 2 (hi). Moreover, we will 
also suppose that F2(x) = Fi(x) = F(x) is an analytic function. These properties generalize 
the known factorization (polynomial factorization if F(x) is a polynomial, see e.g. [2ID 121]) 
properties taking place for the Hermitian case. 

It follows from Q and (@J) that in our case J = o x . 

Keeping in mind the properties of the operators L and J let us introduce the following 
matrix operators: 

*-(2S) MoT)- m 

It follows from the intertwining relations that the operators Qi and Q2 commute with H and 
they apparently are nilpotent. The above factorization properties are equivalent to the following 
anticommutation relation: Q1Q2 + Q2Q1 = F(H). 

Now if we identify our J operator with = r/^ 1 introduced in [201, J = V- — our 
operator with D and JL + J with D% we conclude that the operator Q 2 becomes pseudoadjoint 
to Qx, the operators H,Q\ and Q 2 close a nonlinear superalgebra and one can associate a 
nonlinear pseudosupersymmetry with quantum system described by the Hamiltonian H. In 
the next Section we shall show that a quadratic pseudosupersymmetry may be associated with 
the two-level system. 

4. Intertwining operators for two-level Hamiltonians 

To be able to associate a pseudosupersymmetry with the Hamiltonian given in (j3J) and (j3J) we 
have to find an intertwining operator and a partner Hamiltonian hi. According to Ref. [23] 
the intertwining operator L for a matrix equation such as is defined with the help of a 
matrix- valued function IA =U{t) satisfying the equation 

h U=UL L = diag(/i, l 2 ) (6) 

called the "transformation function" , as follows: 



L = d t - W W = UU~ 



(7) 



Quadratic pseudosupersymmetry in two-level systems 



6 



Here lx and 1% are arbitrary constants. The operator L transforms a solution \l/ of equation 
into a solution $ of the same equation where the matrix Vo is replaced by 

V x = V + AV AV = <yW -W<y . (8) 

Here and in the following the subscript marks quantities before the transformation and 1 
marks these after the transformation. It is not difficult to see that to preserve the form (jlj) of 
the potential so that V\ = ifx^y it is sufficient to take the transformation function of the form 

U=( Ul1 Ul1 V (9) 

^ U 2l -W 2 1 / 

In this case the column-vector Ux = (uxx,u 2 x) T is a solution to the initial equation (jHJ) 
corresponding to the eigenvalue I and the column-vector U 2 = ("Un, — u 21 ) T is a solution to 
the same equation with the eigenvalue —I (note that this symmetry is built into the system 
(J2J)!) so that L in © has the form L = diag(Z, —I). After some simple algebra one finds from 
© that fx = f + Af where 

A/ = l(^-2*)-2/ . (10) 

\u 21 UxxJ 

In general, solutions U 1)2 of equation (j3J) from which the matrix U is composed, U = (Ux, U 2 ), 
are complex, leading to a complex-valued potential difference Af. For physical reasons we 
require real potentials. A necessary condition for Af to be real is that the eigenvalue / be 
purely imaginary. Indeed, it is easy to show that I cannot be real. According to (fTU|) Af is 
defined by the expression — — . Putting — = pexr>(i(p) one finds 

J 1 U21 Mil "21 v ' ' 

Uu M21 t 1 ^ ^ ■< ^ X \ ■ 

= [Q ) COS(y9 + l[Q H ) Sin If (11) 

U 2l Uxx Q Q 

and our claim follows from the fact that g + - is never equal to zero. Finally one can prove 
that I 2 is real (cf. [25]). 

Now when the imaginary character of I is established we see from (II 0J) that the left hand 
side of (fTT|) must be purely imaginary, which is possible only if g = 1, meaning that Uxx and u 2 x 
have the same absolute value. Therefore one can put uxx = pexp(i^i) and u 2 x = pexp(i<£> 2 ). 
Using the fact that Ux = («n, u 2 x) T satisfies equation Q with E = I and setting I = iR, where 
R is real, one gets from (jSJ) a system of equations for p, (fx and (p 2 - Of these equations we need 
only 

02-0! -2/0 + 2^8^(^2-^) =0. (12) 

If R = (|T2|) can be readily integrated. Suppose R ^ 0. The change of the dependent variable 
in equation (|T2j) . ip 2 — <fx = 2arctang, yields for q the Riccati equation 

q + 2Rq-f (l + q 2 ) =0. (13) 

If /o = the equation for q is readily integrated: q = exp(— 2Rt). Considering f Q ^ one can 
linearize (13) by putting q = —u/(ufo), so u is a solution to the second order equation 

u + (2R-f /f )u + f*u = 0. (14) 
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Introducing the new variable if) by putting u = exp(— Rt)y/~foif) one eliminates the first derivative 
term from ()14|) thus obtaining 



ip + 



/o 2 + ^ln/o-[~ln/o 



1 d 
2dt 



R 



ij = o. 



(15) 



This equation has two linearly independent real solutions and, hence, if) is defined up to one 
real constant. Once if) is fixed one calculates q: 



q fo m 

and the potential difference Af = 2i?sin(<£> 2 — (pi) — 2f : 

ARq 



A/ 



l + q 2 



2/ . 



(16) 



(17) 



Solution $ of the equation = E<$> with h x = ~/d t + V u Vx = V + AV, AV = iAfa y 
can be found by applying the transformation operator (JJJ) to solution \I/ of the equation (J3J), 
$ = L^. It is easy to see that the matrix W is diagonal 



W = diag(wi,w 2 ) 



Wi 



-if + Ru 2 i/uu 



w 2 



and the ratio of the components of the spinor JJ\ defining W\ in f)18|) is also expressible in terms 
of the function q: 



u 2 i 

Mil 



(1 + iqf 



(19) 



l + q 2 

Finally skipping calculational details but noticing that just in the same way as it was done 



111 



one can find the following factorizations: 

JL + JL = h 2 -I 2 , LJL + J = h\-l 2 



(20) 



with J = <j\. This means that the function F from Section 3 is F(x) = x 2 — I 2 , the operators H, 
Qi and Q2 close the quadratic superalgebra and the quadratic pseudosupersymmetry underlies 
the two-level system interacting with the electric component of an electromagnetic field. 



5. Application: SUSY transformations of the Rabi oscillations 



In this Section we show a new physical phenomenon we observed while analyzing solutions of 
the system (JTJ obtained using the above developed technique. 

We start with So = 2f = i2 — 0o = constant (this corresponds to the Rabi oscillations (J2J)) 
to get a time-dependent "potential" /i(t) = fo + Af(t) = \ ji\Si{t)t\. Once fi(t) is found we 
calculate the detuning Si (t) = 012 — 0i (t) by integrating the previous equation 

81® = *£h(t)dt. (21) 
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We have found that relatively small but time-dependent perturbations of the field frequency 
0i (t) from its resonance value equal 012 may influence essentially the time behavior of the 
probability -Pi(t) to populate the excited state level with respect to the constant frequency 
case. 

If fo = const equation (fi3j) for ip reduces to 

$ + w 2 tp = w 2 = fl - R 2 = const . (22) 

Solutions of this equation have different properties depending on whether the value w 2 is 
positive, negative or zero. We have found that the oscillating behavior of the probability P\{t) 
disappears when w = 0. In this case the general solution to equation (|22j) is a linear function 
of time ip = At + B which according to gives the following time dependence of the function 
q: q(t) = 1 — A/(Atf + Bf ). Once q(t) is found one calculates the "potential difference" with 
the help of formula (|1U|) and finally the new "potential" / = fi(t): 

m = f °~ 2A*f 2 ti - 2Af (A - 2Bf )t + A* - 2ABf + 2B*f 2 " (23) 

Another restriction leading to the desirable result is A = 2Bf which reduces the previous 
equation to a simpler form 

Since solutions A w (t) and A 2 o(t) of the system ([TJ for / = / = const are known one can find 
solutions A u (t) and A 2 \{t) of the same system with / = /i(t) by applying the transformation 
operator L defined by formulas (|7|). (fTHj) and (fT^j) to the previous solution. In this way imposing 
the initial condition Aq(0) = 1 and ^2i(0) = one finds the probability Pi(t) to populate the 
excited level at the time moment t if at t — only the ground state level is populated 

P 1 (t) = \A 21 (t)\ 2 = 06(1 | 2 4/ft2) [l6/ o 4 0Vcos 2 0t+ 

Afl@t (e - 3f 2 ) sin 20t + (Af%®H 2 + (£ 2 - 3/ 2 ) 2 ) sin 2 0t] . (25) 

Here = a//o + ^ 2 an< i ^0 is the frequency of oscillations of the probability Po(t) (0) at 
/ = fo- It is clearly seen that Pi(t) is an oscillating function provided £ 2 7^ 3/q. For ^ 2 = 3/q 
(0 = 2/o) the probability becomes equal 

p '<" - IT& (26) 

which is a function monotonically growing from zero at the initial time moment till the value 
3/4 at t — > 00. We have to notice that for a fixed £ the parameter f is fixed also, fo = £/\/3, 
which by means of formulas (|24j) and (|2ip fixes the frequency of the electric field in the unique 
way. So, for the given dipole momentum this excitation regime is fixed by the amplitude of 
the electric field. Let us analyze now what is happening with the probability P\{t) when the 
parameters of the model are close to this exceptional point. 
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Suppose now w 2 > and we will consider it to be close to zero. In this case the general 
solution to equation (}2~2~|) may be written as tp = — sinfat+a+b). The function q as given in (JTfij) 
does not depend on the value of the coefficient — but we need this coefficient to realize the limit 
w — > thus recovering the previously obtained solution. Choosing b such that sin 2b = w j fo 
and cos 2b = R/ fo but keeping a arbitrary one gets 

j>_ _ ZJ ZU ~ /o sin(2s7t + 2a) 

This leads to the following expression for q: 

R cos(2zut + 2a) + zu sin(2cc7t + 2a) — f 



fo cos(2zut + 2a) — R 
and finally to the "potential difference" of the form 



(2* 



2w 2 

A ^ ~ Rcos{2wt + 2a) - f ' ^ 

This formula has been previously derived by V.G. Bagrov et. al. by other means [25] ■ Putting 
a = arctg^ — |arctg^ one recovers for = fo + Af(t) the previous result ([23j) as the 

limit — > 0. This means that for w close to zero the probability P\{t) corresponding to the 
potential difference (|2^|) should be close to the previous value The analytic expression for 
P\(t) is rather complicated and we will restrict ourselves by graphical illustrations. 

Let us fix the Rabi frequency 2£. The function Af(t) ()29|1 contains three parameters w, 
fo and a. The parameter f defines the value 20 = 2sJ /q + £ 2 , which is the frequency of 
oscillations of the function Po(t) given by (J2J) to which P\(t) is reduced when the time dependent 
correction A/(£) is absent. As it was already mentioned when fo = (resonance case) the 
function Po(t) oscillates with the Rabi frequency 2£. The parameter w defines the frequency of 
the time dependent correction A/(£) (J29|) for fx = f +Af and the parameter a is responsible for 
the initial value of fi(t). The probability P\(t) is a periodical function if is commensurable 
with w. In this case it exhibits two kinds of oscillations, namely, fast oscillations with the 
frequency 20, which is close to the Rabi frequency when f is close to zero, taking place at the 
background of slow oscillations with the frequency 2w. 

For our numerical illustrations we choose fo — 1. If in standard units this is 1 • 10 11 c _1 this 
corresponds to 10~ n c as the unity of time in our figures. 

Fig. la shows the probability P±(t) for = 2, a = 0.015 and w — 1/4 (solid line) and 
w — 1/6 (dotted line). Fig. lb illustrates the time behavior of the detuning Si(t) calculated 
according to (j2~Ij) for a = 0.015, w = 1/4 and w — 1/6 (solid and dotted lines respectively) 
together with its limiting value corresponding to vj = 10~ 3 and a = 10~ 6 (dashed line). It is 
clearly seen from Fig. la that the period of slow oscillations grows when w decreases and fast 
oscillations go around the limiting value 0.75 with the amplitude increasing with w decreasing. 
Moreover, Fig. lb says that oscillating behavior of P\(t) is transformed into monotonically 
growing one when for w = the detuning becomes a monotone function of time (dotted line 
on Fig lb). If it acquires some oscillating perturbations the probability starts to oscillate also. 
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2.5 5 time 10 12.5 15 -6 



Figure 2a. Probability Pi (t) at Figure 2b. Detuning 5i(t) at 

different values of a. different values of a. 

The next two figures show the dependence of the same quantities on the parameter a which 
is responsible for the phase shift in formula ()29|) at the fixed value w = 1/5. Dotted, dashed 
and solid lines (figure 2a) correspond to a = 0, a = 0.02 and a = 0.08 respectively. Figure 2b 
shows the time dependence of the detuning S(t) for a = (dotted line) and a = 0.08 (solid 
line). From Fig. 2b we can conclude that the parameter a defines mainly the maximum of the 
absolute value of the detuning which it takes at t — 0. Fig. 2a says that the amplitude of fast 
oscillations grows together with a. 




2.5 5 time 10 12.5 15 



Figure 3. Probability Pi(t) at 
different values of 0. 

The next figure shows the dependence of Pi(t) from the frequency of fast oscillations at 
a = and w = 0.2. Dotted, solid and dashed lines correspond to = 2, = 1.6 and = 1.2 
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respectively. More it differs from the critical value equal 2 corresponding to £ 2 = 3/q , when the 
oscillations in formula (|2*5|l disappear, bigger the amplitude of the fast oscillations becomes. 

6. Conclusion 

Using the technique of intertwining operators for a Dirac-like system developed in j2H| we have 
found time dependent electric fields for which the equation of motion for a two-level system 
placed in this field obtained after the rotating wave approximation can be solved exactly. 
Pseudosupersymmetry generators constructed with the help of intertwining operators together 
with the super-Hamiltonian close a quadratic deformation of the superalgebra constructed in 
[2*U] . We conclude, hence, that two- level systems in external electromagnetic fields may have 
hidden quadratic pseudosupersymmetry which is responsible for the new phenomenon consisting 
in disappearance of the Rabi oscillations. 
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